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We consider a pair of adjacent quantum waveguides, in general of different 
widths, coupled laterally by a pair of windows in the common boundary, 
not necessarily of the same length, at a fixed distance. The Hamiltonian 
is the respective Dirichlet Laplacian. We analyze the asymptotic behavior 
of the discrete spectrum as the window distance tends to infinity for the 
generic case, i.e. for eigenvalues of the corresponding one-window problems 
separated from the threshold. 

1 Introduction 

Quantum mechanics exhibits various effects which defy our intuition based on 
"classical" experience. A nice class of examples are bound states in hard-wall 
tubes induced solely by their geometric properties such as bends, protrusions, or 
"windows" . Such systems are interesting not only per se but also from the prac- 
tical point of view as models of various nanophysical devices, and in a reasonable 
approximation also of flat electromagnetic waveguides. 

Among numerous questions such models pose an important one concerns be- 
havior of the spectra in case of two distant perturbations. One can think of it 
as of an analogue of the exponential spectral shift for a pair of distant potential 
wells, despite the fact that the usual methods of the Schrodinger operator theory 
do not work here. The aim of the present paper is to study this problem in a model 
example of a pair of laterally coupled waveguides, or adjacent straight hard- wall 
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strip in the plane, coupled by a pair of "windows" in the common boundary - we 
refer to PP, [2], for a bibliography concerning such models. 

In our recent paper j2] we dealt with the symmetric situation where the widths 
di, d2 of the two channels were the same and so were the window widths ai, 02- The 
technique used in these papers employed substantially the fact that the problem 
can be decomposed into parts with a definite parity, which allows one to study a 
single-window problem with a perturbation which consists of an additional Dirich- 
let or Neumann boundary condition at a segment far from the window. 

The approach based on symmetry works no longer if ai 7^ 02. The man aim of 
the present work is to demonstrate a different technique, suitable for the general 
case, which reduces the question to analysis of a boundary perturbation at the 
distant window. This technique follows the main ideas of where the Dirichlet 
Laplacian in an n-dimensional tube with a pair of distant perturbations described 
by two arbitrary operators was studied. It was assumed in |,4 that these operators 
are defined on functions from W2 vanishing at the boundary, and this assump- 
tion was employed substantially. This is obviously not true in the problem we 
study, since the windows enlarge the domain of the Laplacian beyond the Sobolev 
space At the same time, the general approach of j3] works in our case with 
the appropriate modifications. Moreover, since we restrict ourselves to the two- 
dimensional case and specify the nature of the distant perturbations, we are able 
to obtain a more detailed result in comparison with the general case in j3] . 

In order not to make this study too technical we concentrate in this paper at 
the generic case when the "unperturbed" energy is an isolated eigenvalue of the 
one-window problem, leaving the computationally involved discussion of threshold 
resonances to a sequel. The problem will be properly formulated and the results 
stated in the next section; the rest of the paper is devoted to the proofs. 

2 Statement of the problem and the results 

Let X = (xi,X2) be Cartesian coordinates in the plane, 11+ := {a; : < 0:2 < vr} 
and n~ := {x : —d < X2 < 0}. With the natural scaling properties in mind we 
may suppose without loss of generality that d ^ tt. By 7-1- we denote two intervals 
7± '■= '■ 1^1 =F ^1 < a±! ^2 = 0}, from now on referred to as the windows. The 
numbers a± are assumed to be fixed throughout the paper while the distance 21 
between the windows will be changing playing the role of a large parameter. 

We set n := 11+ U 11" U 7+ U 7_ (cf. Figure 1); the Hilbert space of our 
problem is L2(n). We will employ the symbol H to denote Friedrichs extension of 
the negative Laplacian from the set C^(Il). We will use the symbols cress(-) and 
c"disc(') to indicate the essential and discrete spectrum, respectively. As we have 
indicated in the introduction, this work is devoted to the study of the asymptotic 
behavior of isolated eigenvalues of if as / — +00. In order to formulate the main 
results we have to introduce first some more notations. 

Let fl be an open set in and 7 C fi. Throughout the paper 1^2^(^2,7) 
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Figure 1: Window-coupled waveguides 

will indicate the completion of the set of functions from C°°{fl) having a compact 
support and vanishing in the vicinity of the set 7, taken with respect to the norm 
of the Sobolev space W2{fl). 

We denote 7^ := {x : \xi\ < a, X2 = 0} so that Ua := 11+ U 11^ U 7^ is 
the double waveguide with a single window centered at xi = 0, and Fq := dUa- 
Furthermore, we introduce the corresponding cut-off sets 11^ := Ila fl {x : < b} 
and F^ := F^ fl {x : |xi| < b}. Consider the negative Laplacian in L2(Jla) and 
call H{a) its Friedrichs extension in L2{Ila) from the set C^{Ila) on which it is 
symmetric; by Am(a), m = 1,2,..., we denote the isolated eigenvalues of this 
operator arranged in the ascending order with the multiplicity taken into account. 

The following results were demonstrated in pp, 

Proposition 2.1. For any a > the essential spectrum of H{a) equals [l,+oo) 
while adisc{H{a)) is non-empty consisting of finitely many simple eigenvalues. The 
eigenfunction associated with an eigenvalue Am(a) has a definite parity: it is even 
or odd with respect to Xi if m is odd or even, respectively. In the particular case 
d = TT the eigenfunctions are even in the variable X2. 

The eigenfunctions associated with the eigenvalues A„(a) will be denoted as 
■?/'„(■, a) and assumed to be normalized, i.e. to be unit vectors in L2(Jla)- It is 
easy to check that a) G C°°(Ila). We put a^: := o"disc U a^isc H{a^)] an 
element A* G a^, of this set will be called simple if A.,, belongs to one of the sets 
o"disc(-f^(o±)) only and double otherwise. Furthermore, we set a := (a+,a_). 

With these preliminaries we can formulate the first main result of this paper. 

Theorem 2.1. For any I > 0, a± > the operator H has the essential spectrum 
equal to [1, +00) and finitely many isolated eigenvalues. The number of the isolated 
eigenvalues of H is independent of the window distance provided I ^ max{a_, a+j. 
In the limit I +00 each isolated eigenvalue of the operator H converges to one 
of the numbers from the set cr* or to the threshold of aess{H)- 

By we indicate the set of all bounded domains 5 C having smooth 
boundary and separated from the edges of 7a by a positive distance; we stress that 
the case dS fl dUa 7^ is not excluded. For any A such that Re A ^ 1 we denote 
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Ki = i^tW •= Vl — A and = (A) := — A, where the branch of the root 

is specified by the requirement that the functions are analytic in and y/l = 1. 
The following statements will be proven in Sectional 

Proposition 2.2. In the limit xi — >• ±00 the eigenfunction ipn of H{a) behaves as 

V^„(x,a) = (±l)"+ic(A„,a)e-'^^(^")l^ilsinx2 + 0(e-^"^^l^il), X2 G [0,7r], 
V^„(x,a) = 0(e-R-«r(A)ki|)^ e [-c/,0], 

(2.1) 

if d < TT, and 

^Pnix, a) = (±l)"+ic(A„, a)e-'^^(^")l^'il sin jxsl + O(e-^^^^l^il), X2 G [-vr, vr], 

(2.2) 

m i/ie case of equal-width channels, d = tc. In these relations 

c(A„,a) = ^.^ . / ^„(x,a)e'^i(^")^Mxi = ^ /.^ . / V^„(x, a)e-'^i (""^^^ dxi, 
TTKi (A„) y 7rfi:5^(A„) J 

7a la 

(2.3) 

and c(Ai, a) 7^ 0, i = 1,2. The asymptotic relations \2. i2.^) give rise to valid 
formulcB when both their sides are differentiated. 

Proposition 2.3. For any A G (—00, l)\adisc{H{a)) there exists a unique solution 
of the boundary value problem 

(A + A)f/ = o, xeUa\ja, U = o, xedu^, 



dU 



dx-) 



dU 



X2=+0 8X2 



X2=—0 



belonging to W2(n.a)- For large values of \xi\ this function is infinitely differen- 
tiate and in the limit xi — > +00 it behaves as 



(2.5) 



[/(x,A,a) =c(A,a)e-"^(^")"isinx2 + C(e-^"^^"i), X2 G [0,7r], 

?7(x,A,a) = O(e-^^''^(^)"0, ^2 G [-rf,0], 
if d < Tc, and 

f/(x,A,a) =c(A,a)e-''^(^")^isin|x2| + 0(e-^"^^^^), X2G[-7r,7r], (2.6) 
in the case d = tt, where the coefficient is given by 

c(A, a) = ^7-- [ U(x, A, a)e^^^^^^' dxi. (2.7) 
ttkJ{X) J 

7a 

This coefficient is negative for A < Ai(a). The asymptotic relations \2. ,5j) . 6j) 
give rise to valid formul(E when both their sides are differentiated. 
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For the double window 7+ U 7_ we indicate by S the set of all bounded domains 
C n having smooth boundary and separated from the edges of 7-1- by a positive 
distance; the case dSndU 7^ is again not excluded. For brevity we will introduce 
a two-valued symbol, r := 1 if d < vr and r := 2 if d = vr. 

Continuing the list of the main results we make the following claims. 

Theorem 2.2. Suppose that X^: G a^: is simple being an eigenvalue Xn{a±) of the 
operator H{a±). Then there is a unique eigenvalue of the operator H converging 
to asl ^ +00. This eigenvalue is simple and behaves asymptotically as follows, 

X^{1, a) = X,+ fi^{l, a)e-^'^^(^*)' + C)(/2e-8«^(A.)' + ^-2i{4{x,)+py^^ (2.8) 

^^{l,a) := r7rc(A*, a±)c^(A^,, a=p)/t|(A^,), (2.9) 



where p = p(A) := minj/t^ (A), V4 — A} if d < and p = p(A) := V4 — A if d = t:. 
The associated eigenfunction ip^{x,l,a) satisfies the relation 

/, a) = ^Pnixi T I, X2, a±) + 0(e-2-^(^*)') (2.10) 

in the norms of both the 14^2^ (11) and W^KS") for each 5 G S. 

Theorem 2.3. Suppose that A* G a* is double and A* = A„(a_) = Am(a+)- 
Then there exist either two simple eigenvalues X^{l,a) or one double eigenvalue 
X~{1, a) = X^{1, a) of the operator H converging to X^, as I —>■ +00. The asymptotic 
expansions of these eigenvalues read as follows, 

X^{1, a) = X,± \p{l, a)|e-2«^(^*)' + 0(/e-^"^(^*)' + e'^^^^*)'), (2.11) 
p{l, a) = (-l)'"+V7r/€+(A,)c(A„ a_)c(A„ a+), (2.12) 

Theorem 2.4. Suppose the hypothesis of Theorem \2. !^ holds true. If p{l,a) 7^ 0, 
the eigenvalues A"''(Z, a) X^{l,a) do not coincide and are simple. The associated 
eigenfunctions ilj^{x,l,a) satisfy the relations 

ip^{x,l,a) = ilJn{xi+l,X2,a_)T'ipm{xi-l,X2,a+) sgn/i(/, a) + 0(e"^''^(^*)') (2.13) 

in the norms 0/1^2^^^) and W2{S) for each S* G H. If X~{l,a) = A+(/,a) is a 
double eigenvalue, the associated eigenfunctions ilj^i^x^l^a) satisfy the relations 

^+(x, /, a) = ^„(a;i + /, X2, a_) + ^(e-^''^^"*)'), 
^-(x, /, a) = 7/'„(xi - X2, a_) + 0(e-2«^(^*)'), 

in the norm of 14^2^(11) and W2{S) for each S* G H. Finally, if p{l,a) = and 
X~{l,a) 7^ X^{l,a), the eigenvalues X^{l,a) are simple and the associated eigen- 
functions satisfy the relations 

?/'^(x, /, a) = c^^nixi + I, X2, a_) + c^tp^ixi - I, X2, a+) + Oie''^^^^^*'^^), (2.15) 



+ 



where the vectors := yjtj are nontrivial solutions to the system \4-2(}^ with 
A = A^ such that \\c^\\t3.2 = 1. 
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The leading terms of the asymptotics ()2.8p . ()2.1H) are non-zero provided the 
corresponding coefficients c{X^:,a±) are non-zero. We know from Propositions 12.21 
12.31 that this is true at least for c(A, a) as A ^ Ai(a) or A = A2(a). For instance, if 
Ai(a--) < Ai(a+), the eigenvalue of the operator H converging to Ai(a_) has the 
asymptotic expansion ()2.8|1 . and the coefficient ()2.9|1 of leading term is non-zero. 
Moreover, due to Proposition 12.31 this coefficient is negative. If a± are such that 
Ai(a_) = A2(a+), the eigenvalues of the operator H converging to A* = Ai(a_) = 
A2(a+) have the asymptotics expansions ()2.11|) . and the coefficients of the leading 
terms are non-zero. By Theorem 12.41 the "perturbed" eigenvalues are simple and 
the associated eigenfunctions satisfy the identities ()2.13|1 in this case. We also 
stress that in this case the leading terms of the asymptotic expansions ()2.1H1 have 
the same modulus but different signs. This phenomenon is known in double-well 
problems with symmetric wells. It also occurs in the symmetric case, a_ = a+ and 
(i_ = as we have shown in P|. 

We conjecture that the coefficient c(A, a) is non-zero for all values a and A < 1. 
If it is true, this fact would imply that the leading terms in the asymptotics ()2.8j) . 
()2.11|) are non-zero. In turn, this fact together with Theorem 12.41 would imply that 
a double A,,, G a* splits into two simple "perturbed" eigenvalues and the formute 
()2.13|) are valid for the associated eigenfunctions. 



3 Analysis of the one-window problem 

In this section we shall study the following boundary value problem, 
(A + A)n = 0, X ella\ja, u = 0, X e Ta, 



du 



dXri 



du 



dXn 



/, X e -fa- 



(3.1) 



The function / is assumed to be an element of L2(7a)- A solution to this problem 
is understood in a generalized sense, more specifically, as a function belonging to 
iy2Hna, r^) for each 6 > and satisfying the equation 

-{Vu, VC)L,(n.) + A(n, C)L,(n„) - (/, C)l,(,„) = (3.2) 

for any ( G (!!„). By standard smoothness- improving results about solutions 
to elliptic boundary value problems, cf. [S"; Ch. 4, §2], the said solution belongs to 
C°°(n+ U 11" \ 7^). As we have said in the introduction we will deal in this paper 
with the non-threshold case only. Thus the parameter A is supposed to belong to 

for a fixed 6 > 0, where is a set of all A separated from the halffine [1, +cxd) 
by a distance not less than S. 

We seek a solution to the problem ()3.1|1 belonging to L2 (!!„). We fix jS > a 
and put P := {x : \xi\ < a + (3, < X2 < do}. The number do here is chosen so 
that do < d, and the lowest eigenvalue of the negative Laplacian in P subject to 
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Dirichlet boundary condition on dP \ 7^ and to Neumann one on 7^ exceeds two. 
We consider the boundary value problem 

du 1 

(A + A)^I = 0, xeP, u = 0, xedP, ^ = 2^' ^2-2) 
which is again treated in the weak sense, 

-{Vu, VC)l.(p) + X{u, Omp) - \{f, Ol,m = (3.4) 

for each function ( G C'^{P) vanishing in a neighborhood of 9P\7q. The problem 
(jH.Hj) is uniquely solvable in the space Vl^2^(P, dP \ ■ja) and the solution belongs to 
C°°(P \7j - see [3 Chap. E, §5, Rem. 5.1] and Chap. IV, §2]). 

Let Xi = Xi{^) be an infinitely different iable function, even w.r.t. the variable 
X2, equal to one if < a + /?/6 and \x2\ < do/ 6, and vanishing for > a + P/3 
or 1x21 > do/3. We extend the function u in an even way for 0:2 < setting 
u{x) := u{xi, —X2) as 2:2 < and denote Uf{x) := Xi{^)u{x). 

Lemma 3.1. The function Uf belongs to W2(Jla,^a) H C°°(Ila) and satisfies the 
equation 

-(Vuf, VC)L2(n„) + A(m/, C)L2(n„) - (/, Oi^i^a) = OiaCn.) (3.5) 

for any ( G C^(Jla), where 

F = Ti(A, a)f := 2Vu ■ Vxi + uAxi . 

The operator Ti : -^2(70) L2{{x : |xi| < a + P/3,\x2\ < do/3}) is linear, 
bounded, and holomorphic in X. The operator T2{\, a) f := Uf is linear, bounded, 
and holomorphic in X as a map from ^2(70) into W2(Jla, ^a), ^^2^(5*), and VF2^(n^\ 
nf), where S is such that S C 11^ or S C . 

Proof Let {u^^^ be a sequence of functions from C°°{P) vanishing in a neighbor- 
hood of (9P\7„, which converges to u in 14^2^(P). It is easy to see that the functions 
u^j\x) := xi{x)u^^\x) belong to 14^2^(11(1, F^), and that they converge to Wf in the 
norm of 1^2^ (Ila) as j cx3, so G 1^2^(110, F^). Next we observe that Uf belongs 
to C°°(Ila) as it follows from the fact that u G C°°{P \ 7,^). Since the function 
Uf is even in the variable X2, we find that for each ( G C^{Ila) the left-hand side 
of ()3.5|) equals twice the expression 

-{VUf, VC+)l.(P) + X{Uf, C+)l,(P) - (/, C)L,i,^) , 

where C^{x) := ({x) + ({xi, — X2). In view of ()3.4|1 and the definition of Xi "we get 
- (Vm/, VC)L2(n„) + X{uf, C)L2(n„) - (/, C)L2(7a) 
=2 1^ - {Vu, V(xiC+))l2(p) + X{uf, XiC)L,iP) - lif, XiC+)l,(>) 
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+ (Vn,C+Vxi)L,(P) - (^IVxi, VC+)l.(p) 

=2((Vm,C+Vxi)l,{p) - (wVxi, VC+)l,(p)) 
=2((Vm,C+Vxi)l2(p) + (divMVxi,C"')L2(P)) = {F,OL,in^)- 

The boundedness of the operator Ti(A, a) follows from the above mentioned theo- 
rems on improving smoothness of solutions to elliptic boundary value problems. 

In order to check that Ti is holomorphic in the variable A we just need to show 
that the mapping / i-^ m is bounded and holomorphic as an operator family from 
L2(7a) into W^{P) and W^{S f] P), where ^ G S„ and ^ C TF. To prove the last 
claim it is sufficient to reduce the boundary value problem to an operator equation 
in W2{P, dP \ in the standard way - see Ch. I, §2] - and to apply then 
Proposition 4.5 of ^, Ch. XI, §4]. □ 

We seek the solution to the problem ()3.1|) in the form u = uj + u. As is 
compactly supported, the function u has to be an element of L2(Jla)- It follows 
from ()H.2j) . that the function u must also obey the integral relation 

-(Vn, VC)L.(n.) + A(n, C)L.(n.) = -{F, C)L.(n.) (3.6) 

for any ( G C^(Jla). Thus u has to solve the boundary value problem 

-(A + A)m = F, xeUa, u = 0, xeTa, (3.7) 

belonging to L2{Ua) and W^{U''^,rl) for each b > 0. By Theorem 4.6.8 of [HI 
Ch. 4, §4.6] any solution of this problem belonging to L2{Ila) is an element of the 
operator domain of H{a). In this way the problem ()3.7|) can be cast into the form 
{H{a) — X)u = F, which in turn gives u = [H{a) — A) ^F. 

Let us next denote T3(A,a) := T2(A,a) + {H{a) — A)^^Ti(A,a). In order to 
analyze properties of this operator we need an additional notation and a lemma. 
For any numbers &i, 62,^3 G M we set fl± := {x : ±Xi > bi, 62 < 3^2 < ^3} and 
LU± := {x : ±Xi > 61} n dil±. 

Lemma 3.2. Let v G W2{^±) be a solution to the problem 

(A + A)w = 0, xe^^, v = 0, xeuj^, 
and < 62 — &3 ^ TT, A G S^. Then the function v can be represented as 



f (X, 
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where 

63 

2 [ nj 

(^iW '■= 1 / v{ai,X2, X) sin- — (x2-&2)dx2. (3.9) 

63 - 62 J 03 - 62 

The series \'^. tSj) converges in the norms ofW^{{x : ±Xi > 64, &2 < X2 < ^3}), 
m ^ 0, /or any 64 > 61. T/ie coefficients aj satisfy the condition 

00 

2 XI = ll^(±^i> ■> ^)IU2(b2,f.3)- (3-10) 

This lemma is a particular case of Lemma 3.3 of 4J so we skip the proof. 

Lemma 3.3. The operator T3(A,a) is bounded and meromorphic in X & Ss as 
a map from -^2(70) into W2(Jla,^a) and into W2{S) for each S E Ea- Its poles 
coincide with the eigenvalues of the operator H{a) . For any X close to an eigenvalue 
Xn of H{a) the representation 

T3(A, a) = Y^^4(a) + T,{X, a) (3.11) 

A — An 

holds true. Here Ti{a)f := (/, V'n)-L2(7a) ^.nd the operator T^ is hounded and holo- 
morphic in X E E>s as a map from -^2(70) into W2(Jla, ^a)- The operator T5 is also 
bounded and holomorphic as a map into W2{S) for each S E Sa- 

Proof. In accordance with jHl Ch. 5, §3.5] the operator {H{a) — X)~^ is bounded 
and meromorphic in L2{Ila), its poles coincide with the eigenvalues of H{a) and 
for A close to A^ the representation 

{Hia) -X)-' = -^-i!^(.,^J^^(n.) +T6(A,a) (3.12) 

is valid, where the operator Tq{X, a) is bounded and holomorphic in A in the vicinity 
of A„. The function u := TQ{X,a)F is a solution to the boundary value problem 
()3.7p with F replaced by F := F — {F, ipn) L^(^iia)'^n] it means that 

l|V^lli2(n„) - A||M||L(n„) = {F,iL)L,(Tia)- 

This relation together with ()3.12|) imply that the operator Tg is bounded and 
holomorphic as a map into 1^2^(110) as well. Using again the smoothness-improving 
theorems mentioned above we conclude that the operator Tg is also bounded and 
holomorphic in A as a map into W2{S) for each S G H^. 

Since the function ipn is an element of W2(Jla,^a), the relation ()3.5|1 is valid 
for ( = ipn- For any / G L2{'ja) the function Ti(A,a)/ is compactly supported, 
hence we have 

(Ti(A„,a)/,V'n)L2(n„) = -(Vu/, V?/'„)L2(n„) + A(n/, ?/^„)L2(n„) - if,^n)L2M- 
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According to Lemma EUl the function Uf belongs to W2(Jla,^a), which allows us 
to proceed with the calculations, 

(Ti(A„,a)/, V'n)L2(n,) = -(/, ^n)L2(7a) = -Tif- 

Substituting the relation thus obtained together with ()3.12|) into the definition 
of the operator T3 and taking into account Lemma we arrive finally at the 
statement of the lemma. □ 

Let us next fix a number a > 0. For any / ^ (a + a) we define operators 
T^(A, /, a, a) which map an arbitrary v G VF2^(nj^) into the function 



(A)(xiTa)g-24(A)« 



^^±gT«:-(A)(xiqia)^-2K-(A)Z 



e J 



a 



v{a,X2)smjx2dx2, jSf 



via, X2) sin -^^2 dx2 



-A, J ^2. 



(3.13) 



The branch of the root in the definition of the functions kj is specified by the 
requirement that the functions are analytic in and a/T = 1. 

Lemma 3.4. The operators : VF2^(n") -^2(75) o'^e well defined, hounded and 
holomorphic m A G §5. The estimates 



d'T^ 



^ Cre~(2«-a-a)ReK,{A)^ i = 0, 1, 2, 



hold true uniformly w.r.t. A G S5 and I ^ {a + a). 

Proof. We will prove the lemma for T^ only, the argument for T^ is similar. The 
function u belongs to 14^2^(111^), hence we have the estimate 

00 

where the constant C is independent of A G and / ^ (a + a). Employing this 
inequality we infer that 



^2(75) 



~{2l-a~a)ReKj{X) 



,=1 JRe4(A) 
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0=1 



^''^ ^ °° A2^-2{2l-a-a)ReK+{X)\ 



0=1 



l<(A)| 



j2 exp (-2(2/ - a - S) Re(K+(A) - (A))) ' 
v/|ImA| 



1/2 



where C is independent of A G S^. In the same way one can prove that 



(A)(-a)g-2K-(A)Z 



^ (2i— a— a) Re (A) 1 



L2iTa) 



The last two estimates imply that the operator T^" : 14^2^(11^) -^^2(72) is well 
defined and bounded. One can check easily that 



2/€T(A)d 



Repeating the argument which yielded the estimate for T^v we can establish that 

^ (2i— a— a) ReKr(A) ||„,|| 



dT:^v 



dX 



L2{Ta) 



with the constant C independent of A G §5 and I ^ (a + a). Consequently, the 
operator exists, it is bounded and the stated estimate for its norm holds true. 
The norm estimate for is obtained in a similar way. □ 



For any / ^ (a + a) we define operators T^{\^ I, a, a) which map any / G L2{'ja) 
into the function 

On dzL 

'xi ± 2/, +0, A) - 7— (xi ± 21, -0, A), xi e {-a, a). 



dx 



dx2 



Here m is a solution to the boundary value problem ()3.ip belonging to M/2^(n). 
Taking into account Lemma f3. 31 together with the boundedness of the embedding 
PF2^(na) into L2{{x : \xi ± 2l\ < a,X2 = 0}), we conclude that the operators 
: Iv2(7a) -^2(72) are bounded and holomorphic in A G §5. 

Lemma 3.5. The poles of the operators coincide with the eigenvalues of the 
operator H{a). For any compact set K C Ss separated from crdisc(-f^(«)) by a 
positive distance the estimates 



< Cre-2'^"'^^(^), i = 0,l, 



(3.14) 
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hold true with C which is independent of X & K and I. For any X close to an 
eigenvalue A„ of the operator H{a) the representation 



Tg (A, /, a, a) 



A -A, 



-T^i^a) + T^{X, I, a, a) 



(3.15) 



is valid, where 



dtp dijj 
0^(a;i, /, a) := — -^(a;i ± 2/, +0, a) - — -^(xi ± 2/, -0, a), xi G (-a,a). (3.16) 
0x2 0X2 

The operators Tg : Iv2(7a) — * -^2(72) o-f& hounded and holomorphic w.r.t. X in the 
vicinity of A„ and satisfy the estimates 



dX' 



^ ^;^+lg-2me.+ (A)^ ^ = 0,1, 



(3.17) 



where the constant C is independent of X and I. 
Proof. Due to Lemma [3.21 we have 

T^(A, I, a,a)f = T^(A, /, a,a)u. 



(3.18) 



Here m is a solution to the boundary value problem ()3.H1 . Using this identity 
and the representation ()3.11|1 . we arrive at ()3.15j) . where Tg"*" is bounded operator 
holomorphic in A. Moreover, 



T^(X, I, a, a) 



T^{X, I, a, a) — T^[Xn^ I, a, a) 



X — Xr, 

X 



T4(a) + T.^{X, I, a, a)T^{X, a) 



1 



dX 



(A, Z, a, a) 



X — Xn J dX 

An 

A A 

1 



z, I, a, a) dz T4(a) + T^{X, I, a, a)T^{X, a), 



(A - XnY 




dX^ 



Z2,l,a,a) dz2dzi T4(a) 



A„ zi 



d 



+ — {T,^{XJ,a,a)T^{X,a)). 



Applying now Lemma [3.41 we obtain the estimates ()3.17j) . 

The operators , i = 0, 1, are bounded uniformly in A G i^, thus in view of 
the relation ()3.18j) and Lemma [3.41 we arrive readily at the estimates ()3.14j) . □ 

Concluding this section we shall prove Propositions 12.21 and 12.31 

Proof of Proposition \2/A Applying Lemma [3. 21 to ipn with bi = ±a, 62 = 0, 63 = vr 
and 62 = -d, 63 = 0, we obtain the formula; ^i^, (Q- The factor (±1)"+^ in 
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these formulae is due to the definite parity of ipn w.r.t. Xi. The formula ()2.3p for 
c["^ follows from the chain of relations obtained by integration by parts, 

0^/e-f(-.-.si,>..(A + A„),„(..a)d. 

n+ 

' e±''^(^")"^V'n(x, a) dx - (±1)"+Vk+(A„)c(A„, a). 

It remains to check the inequalities c(Aj,a) 7^ 0, z = 1,2. The eigenfunction ijji 
associated with the ground state can be chosen non-negative. Moreover, ipi is not 
identically zero at 7a, since otherwise it would be an eigenfunction of the negative 
Dirichlet Laplacian in n"*" and would correspond to the eigenvalue Ai < 1. At the 
same time, the spectrum of the mentioned operator is the halfiine [l,+oo). The 
described properties of ipi and the formula ()2.12|1 imply that c(Ai, a) 7^ 0. 

According to Proposition 12.11 the eigenfunction ip2 is odd w.r.t. Xi. It allows 
us to modify the formula ()2.3j) . 

c(A2,a) = — ,\ , / ij2ix, a) smh fit {X2) X I dx I 
TTKi (A2) J 

2 



7a 

a 



2 



(3.19) 



4'2{xi, 0, a) sinh Kf{\2)xi dxi. 







The eigenvalue A2 is the ground state of the negative Laplacian in Ilar\{x : Xi > 0} 
subject to Dirichlet boundary condition on dUa (1 {x : xi > 0} and to Neumann 
one on {x : Xi = 0, —d < X2 < vr}. Hence 4'2{xi, 0, a) ^ 0, and 4'2{xi, 0, a) ^ as 
xi e (0,a), and by ()3.19j) these inequalities imply that c(A2,a) 7^ 0. □ 

Proof of Proposition \2.'A The unique solvability of the problem ()2.4j) is ensured 
by Lemma f3. 31 Moreover, we have U = T^{\,a)e'^"^'^^^^'^ . The relations ()2.5|) . ()2.6p 
follow from Lemma f3.2| and the formula ()2.7|1 is proved in the same way as ()2.3j) . 

Integrating by parts and employing the formula ()2.7j) . we obtain a chain or 
identities, 



= y f/(A + \)Udx = A||f/||2^(n) - ||Vf/| 



2 



U 



dU 



_ dU_ 

X2=+Q 0X2 



dx 

la 

=A||f/||L(n)-||Vf/||i,(n)-vr«:+(A)c(A,a) 



X2 = -0 



dxi 



which implies 



c(A..)^^"^"^""'Lr"^''"'. (3.20) 
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Since U G (11, 911) , the minimax principle yields the inequality 
We substitute this inequahty into the formula ()H.2()|1 and obtain 

^(^■"'«^(a^-0 liv;/|iL,n, < 0, 

ifA<Ai(a). □ 

4 Reduction of the perturbed problem 

After this preliminary let us turn to our main problem; we are going to reformulate 
it as a suitable operator equation. Recall that we are looking for eigenvalues of 
the operator H , i.e. non-trivial L2(n)-solutions to the boundary value problem 

-Alp = Xip, xeu, ip = o, X e an. (4.1) 

We denote := {x : —b < xi < b, —d < X2 < n} and introduce the cut-off regions 
n'' := := annQ^ Solutions to the problem (jUD) can be identified with 

functions belonging to 1^2^(11^, T^) for any 6 > such that 

(VV^,VC)L.(n) = A(V^,C)L.(n), (4.2) 

holds for each ( G C^(n); it follows from the smoothness-improving theorem 
mentioned above that such a ip belongs to C°°(n). 

We assume that A G S^, with 5 > is chosen in such a way that cr^, C S^. 
Let f± = f±{-,l) G L2(7a±) be an arbitrary pair of functions. Denote by u± the 
solutions of the problem 1)3.11) with a = a± and f = f± & L2(7a±) and assume that 
u± G L2(Jla±)- We will seek a solution to the problem ()4.1|) in the form 

ip{x, A, /) = M+(xi — /, X2, A, /) + ■u„(xi + /, X2, A, /). (4.3) 

Suppose for a moment that the function defined in this way solves the problem 
()4.H) . In such a case the function ijj is infinitely differentiable at the points of the 
segments 7^^, and therefore 

dtp dip 

— (xi, +0, A, /) - ^(^b -0, A, /) = 0, X e 7±. 

Substituting from ()4.3|) into this identity, we obtain a pair of equations, 

f±{xi) + -^{xi±2l,+0,\J)--^ixi±2l,-0,X,l) = 0, xG7a±. (4.4) 

Denote / = (/+,/„) G i^2(7a+) © -Z^2(7a_)- The following lemma states that the 
last equation is equivalent to the original problem ()4.1)1 . 
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Lemma 4.1. To any solution f G i^2(7a+) © -^2(7a_) of ^-41 ) O'nd functions u± 
solving for a = a±, f = f± there exists a unique L2(J1) -solution of given 
by \4-^ - Reversely, for any solution tp of \4-l^ there are unique f G 1^2(70+) © 
-^2(7a_) solving \4-4\ ) o^^'^d unique functions u± G L2(Jla±) satisfying \S. 1\) with 
O' = ci±, f = f± such that ip is given by \4-3}) - This equivalence holds for any 
A G §5 and I ^ max{a_, a+} + 1. 

Proof. Suppose that / G i^2(7a+) © -^2(7a_) is a solution to the equations ()4.4|) . 
where the functions u± G L2{Ila±) solve the problem p.ip a = a± and f = f±- We 
define ip in accordance with ()4.H|1 . The functions u-t are elements of ^2(1!), hence 
the same is true for ip. Moreover, the function belongs obviously to 1^2^(11'', T^) 
for each b > and vanishes on F. 

Let us check that the function tp satisfies the equation ()4.2j) . To this purpose, 
we indicate by X2 = X2{xi) an infinitely differentiable cut-off function being equal 
to one if |xi + /| < max{a+, a^} + 1/2 and vanishing if \xi + /| > max{a_|_, a_} + 1. 
For any ( G (7^(11) we have 

(Vm+(xi - /,X2, A,/), VC)^^(n) ~ ^{u+{xi - I, X2, X,l), C) L^(^n) 
= {Vu+{xi - /,a;2, A,/), V(Cx2))^^(n) - ^{u+{xi - l,X2, XJ)XX2)) ^^^^^^ 
+ {Vu+{xi - /,a;2, A,0, V(C(1 - X2))) ^^^^^^ - A(m+(xi - /, X2, A, /), C(l - X2)) ^^^^^y 

(4.5) 

Since u+{-) is an element of C°°(na^ \Q +), we can integrate by parts, 

(Vm+(xi - /,X2, A,/), VCx2))^2(n) - A(m+(xi - l,X2,X,l),CX2) ^^^^^^ 
^(xi - /,+0, A,/) - ^{xi -0,A,/) Cdxi- 



(7X2 c'a;2 
- / CX2(A + A)M+(a;i - Z,X2, A,/) dx = / + /)Cdxi. 



n 



We have employed here the equation satisfied by u+ as well as the equation ()4.4j] 
for /+. Since C,{xi + l,X2){l — X2{.Xi + l)) G C^ijla+)i we can use the identity ()3.2p 
to infer that 

(V?i+(xi -/,X2, A,/), VC(1 - X2))i2(n)-(«+(a;i - ^ a;2. A, /), C(l - X2))^^(n) 

= - / f+{.xi - /)Cdxi. 



7+ 



We substitute now the last two relations into ()4.5|1 and arrive at the identity 
(Vm+(xi - /,X2, A,/), VC)^^(n) ~ ^("+(^1 - ^'^2, A,/),C)^2(n) 

= U-{Xl + /), C)l2(7-) - U+{Xl - 0' C)l2(7+)- 
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In the same way one can check that 

(Vu_(xi + /,a;2,A,/), VC)^^(n) ~ + l,X2, X,l)X) ^^(^yi) 

= if+i^i - 0> C)l2(7+) - if-{xi + I), C)l2(7-); 

summing the last two relations we arrive at the relation ()4.2p for the function 

Let be a solution to the problem ()4.H1 belonging to L2(n). By smoothness- 
improving theorems the function ip belongs to C°°{{x : —1 ^ Xi ^ 1,0^ X2 ^tt}) 
and to C°°{{x : — 1 ^ Xi ^ 1, —d ^ X2 ^ 0}). This allows us to define the numbers 

= "^"(-^5 := - y 1^^(05 X2, A, /) ± ^^(0' 2:2, A, l)j sin jx2 da;2, 

(5f = /3f (A, ^)--=\] (^^(0' ^2, A, /) ± ^^(0, X2, A, /) j sin "^x^ dx2. 

Using these numbers, we introduce the functions u± in the following way: 

00 

=F ^ a^2, A, /) := a^(A, /)e^'^j sin JX2, ±Xi^O, X2G(0,7r), 



00 



m±(xiT/,X2,A,/) := ^/3+(A,0e^''^ ''^sin^a;2, ± Xi ^ 0, X2e(-rf,0), 

u±{xi =F ^ 3^2, A, /) := ^/'(x, A, /) — u±{xi ± /, X2, A), ± Xi > 0, X2 G (— d, vr). 

Proceeding in the same way as in the proof of Lemma 4.1 in ^j, we check that the 
functions u± are well defined and 

u±eWl{Tl,^Ja^)^Wi{S), SeEa, (4.6) 
(A + A)m±(xi =f/,X2,A,/) = 0, X G n\ {x : xi = 0}. (4.7) 

The relation ()4.3|) follows from the definition of the functions u±. Now we set 

/±(xi, /) := -^{xi ± 21, +0, A, + -^(xi ± 21, -0, A, I), Xi G (-a±, a±); 

CX2 CX2 

(4.8) 

in view of ()4.6|1 we can conclude that f± G L2 (7a± ) • We also note that the definition 
of u± and the smoothness of ip at 7-1- imply 



/±(xi,/) = — ^(xi,+0,A,/) - — ^(xi,-0,A,/), xiG(-a±,a±). (4.9) 

CX2 CX2 



Let us check the integral equation ()3.2|) for the function = m+(x, A). Taking 
into account ()4.7|) . ()4.9|) . and integrating by parts, we get 

- (Vm+, VC)L2(n„) + A(m+, C)L2(n,) 
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+0, A) - -0, A), c) = (/+, C)l.(,.,). 

for any G C^{U.a+)- In the same way one can check that 

-(Vm_, VC)L2(n,) + A(M_,C)L2(n,) = (/-,C)L2(7a-) 

for any ^ G C^(Jla_), thus m± are solutions to the problem ()3.1|) with a = a± and 
f = f±. The equations (jO)) follow from (jOI)- □ 

Suppose that \ E E>s\a^:. In that case the functions u± introduced above can 
be represented as u± = T3{X,a±)f±, thus the equations ()4.4|) become 

/ + T8(A,/,a)/ = 0, (4.10) 

where the operator Tg : i^2(7a+) © -^^2(7a_) ^ -^2(7a+) © -^^2(7a_) is defined by 

Ts{X,l,a)f := (r8+(A,/,a_,a+)/_,Tg-(A,/,a+,a_)/+). 

Now we are ready to demonstrate the first one of our main results. 

Proof of Theorem \2.1\ If a± = the essential spectrum of the operator H is obvi- 
ously [l,+oo), and an elementary argument using Dirichlet- Neumann bracketing 
[Tin Ch. Xni, §15] and the minimax principle piIl Ch. XEI, §1] shows that the 
threshold of the essential spectrum of H is one, i.e. a^ss^H) C [1, +oo). The oppo- 
site inclusion can be shown easily; one needs to employ Weyl's criterion (see, for 
instance, proof of Lemma 2.1 in j^). 

The operator H being self-adjoint, its isolated eigenvalues are real, and in view 
of the above observation they are smaller than one; we arrange them conventionally 
in the ascending order counting multiplicity. Next we use bracketing again in a way 
analogous to P : we add Neumann boundaries at segments corresponding to Xi at 
the endpoints of 7± and X2 G (— (i, vr). In this way we get an operator estimating 
H from below, and since only the window parts contribute to the spectrum below 
one we infer by minimax that H has finitely many eigenvalues for any / > and 
their number has a bound independent of /. 

Let K C §5 be any compact set separated from cr^, by a positive distance. By 
the estimates ()3.14|) the operator Tg has a norm being strictly less than one for 
A G X and I large enough. For such A and I the equation ()4.1()|1 has thus a trivial 
solution only, and in view of Lemma 14.11 this implies that the operator H has no 
eigenvalues in the set i^' if Hs large enough. This means that each eigenvalue of 
the operator H has to converge to one of the numbers from the set a^, or to the 
threshold of the essential spectrum. □ 

The eigenvalues if, i.e. those A for which the problem ()4.H1 has a nontrivial 
L2(n)-solution, coincide in view of Lemma (4.11 with the values of A for which the 
equation ()4.4|) has a nontrivial solution. In the case considered here we deal only 
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with the eigenvalues of H which converge to a value A,,, G a* separated from the 
threshold, in other words, being smaller than one. 

Our aim is to solve the equation ()4.4|) and to obtain in this way an equation 
for the aforementioned values of A. Consider a A* G 0"=,,; if A=k = A„ is an eigenvalue 
of the operator H{a+) we set 

(Pti-J) ■■= (O,0;(-,/,a+)) G L2(7aJ ©i^2(7a_), T+/ := (/+, 

where 0~ is determined by ipn in accordance with ()3.16|) and ipn is an eigenfunction 
associated with A„, in the opposite case we put 

:= (0,0) G L2(7a+) ©^2(7a.), T+f := 0. 

Analogously, if A* = A„ is an eigenvalue of H{a^) we set 

0;(-,O:= (0+(-,/,a_),O) GL2(7aJ©i:2(7a_), T^"/ :=(/_, 

where (j)^ corresponds to ipn according to ()3.16|) and ipn is an eigenfunction asso- 
ciated with Aji, otherwise 

:= (0,0) G L2(7a+) ©i^2(7a_), T^f := 0. 

Given a number A* G cr*, we consider the equation ()4.4|1 for A in the vicinity 
of A=K. Assume first that A 7^ A,,,, in which case the equation ()4.4j) is equivalent to 
()4.10p . In view of Lemma 13.51 the operator Tg is bounded and meromorphic as a 
function of A G §5, and the numbers A,,, G 0"=,, are poles of Tg. For any A close to A^, 
the operator Tg can be thus represented as 

T8(A, /, a) = cf>t{; l)-^ + C(-, + T9(A, /, a), (4.11) 

A — A* A — A* 

where the operator Tg acts as 

Tg{X,l,a)f := (Tg+(A,/,a_,a+)/_,T9"(A,/,a+,a_)/+) 
if A* G (Tdisc(-f^(a+)) \ adisc(-f^(a-)), 

Tg{X,l,a)f := (T9+(A,/,a_,a+)/_,Tg"(A,/,a+,a_)/+) 
if A* G (Tdisc(-f^(a-)) \ o-disc(-f^(a+)), and finally, 

Tg{X,l,a)f := {T^{X,l,a_,a+)f^,Tg{X,l,a+,a^)U) 

if A.^ G adisc{H{a_)) n adisc{H{a+)). The operator Tg on T2(7a+) © T2(7a_) is 
bounded and holomorphic w.r.t. A in the vicinity of A,,,, and the estimate 

^ Cr+^-^^^^^tW^ i = 0,l, (4.12) 



dX^ 
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holds true with a constant C which is independent on A and /. 
We substitute the representation ()4.11|) into ()4.10p to obtain 

A — A — A* 

Since the norm of Tg is small for large / due to ()4.12p . the operator (I + Tg)"^ is 
well defined being bounded in L2(7a+) © -Z^2(7a_)- We apply this operator to the 
last equation arriving at 

A — A* A — A=K 

where $f (■, A, /) = (! + Tg(A, /, a))^^(pf{-, I). The last equation implies that 

/ = c+$+ + c_$; (4.14) 
for some numbers c±. We substitute from here into ()4. 1H|) obtaining 

(4.15) 

where the quantities Aij = Aij{X,l) are defined by 

AniX, I) := T+$+(-, A, /), A2(A, /) := T+$;(-, A, /), 
A2i{\, I) := T^Ki; A, /), A22{\, I) := r4-<l>;(-, A, /). 

The definition of together with the estimate ()4.12j) imply for / large enough 

*f = 0f + O(e-''^^''^(')il0f||). (4.16) 
If </)^ 7^ 0, and = 0, in particular, we have 

$+^0, $;=0, ^12 = ^22 = 0, (4.17) 
and in this case the equation ()4.15|1 holds if and only if 

If / corresponds to an eigenfunction ip of the problem ()4.ip by ()4.3p . the number 
c+ is non-zero. Indeed, in the opposite case ()4.14p and ()4.17|) would imply that 
/ = 0, which by Lemma Wl\ results in ip = 0. Consequently, the equation ()4.1Up 
has in this nontrivial solution if and only if 

A-A, + An(A,/) =0. (4.18) 



19 



If A is a root of this equation, the corresponding nontrivial solution of ()4.10|) can 
be expressed as ()4.14|) with c+ 7^ and c_ = 0. 

In the case 0^" = and (f>~ 7^ similar arguments lead us to the conclusion 
that the equation ()4.10j) has a nontrivial solution if and only if 

A-A, + A22(A,/) =0, (4.19) 

and the corresponding non-trivial solution can be written as fj4.14|) with the coef- 
ficients c+ = and c_ 7^ 0. 

Finally, if both the functions (f)f are non-zero, they are linearly independent 
by definition and the same is true for the functions $f . Hence the equation ()4.1()j) 
holds if and only if 

((A-A,)E + A(A,/))c = 0, (4.20) 
where E is the unit matrix, and 

AM n - AuiXJ)\ ^._(c, 

^22(A,/)J' 

The column c is non-zero, since otherwise ()4.14j) and ()4.17j) would imply / = 0, 
thus the system ()4.20|) of linear equations has a nontrivial solution if and only if 

det((A- A,)E + A(A,/)) =0, (4.21) 

which can be rewritten as 

(A-A,)2 + (A-A,)trA(A,/) + detA(A,/) = 0; (4.22) 



the corresponding non-trivial solution of the equation ()4.10|1 is given by ()4.14|1 . 
where {Vl) is a nontrivial solution of ()4.2()|1 . 

Assume now that A = A*. Let A* coincide with an eigenvalue A„ of the operator 
H{a^) being not at the same time an eigenvalue of H{aJ). In this case we again 
can claim that m_ = T3(A=|,, a_)f_, on the other hand, the boundary value problem 
for with A = A* is solvable if and only if 

0= J /+^„(a;,a+)dx = T+/. (4.23) 

This follows from Lemma IIOI The function m+ is given by u+ = T5(A=k, a+)/+ — 
c+ '?/'+, where c_|_ is a constant. We can substitute now the described u± into ()4.4j] 
and obtain 

/ + T9(A,/,a)/ = c+0+, 

/ = c+$+. (4.24) 

This function will generate a solution to the problem ()4.H1 if and only if ()4.18j] 
holds true. Substituting ()4.24j) into ()4.23|1 . we arrive at the equation ()4.18j) with 
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A = A*. If c_|_ 7^ holds in ()4.24|) we see that the formula ()4.23|) coincides with 
()4.14|) with c_ = 0. Consequently, in the case A* G a{H{aJ)) \ a{H{a^)) the 
equation ()4.18|) determines all the values of A in the vicinity of A* for which the 
equation ()4.4|1 has a nontrivial solution; these nontrivial solutions are given by 
with c+ 7^ and c_ = 0. 

In the same way one can check that the equation ()4.19j) determines the sought 
values of A in the case when A.,, is an eigenvalue of the operator H{a_) and not of 
H{a^). The corresponding nontrivial solutions of ()4.4|) have c_|_ = and c_ 7^ 0. 

Finally, if A* G cr* is double and A = A*, the solvability conditions of the 
boundary value problems for u are T^f = 0. If this holds true, the functions u± 
are given by u± = T5(A,a±) — c±ip±{-), where c± are constants and iIj± ^t^^ the 
eigenf unctions of H{a±) associated with A*. The equation ()4.4|) becomes 

/ + Tg{X^, I, a)f = c+(j)t + c_0;, 

which yields the relation ()4.14j) . The solvability conditions f = are nothing 
else than the system of linear equations ()4.20j) . In this way ()4.14|) . ()4.20|) . and 
1)4.2111 describe the sought values of A in the vicinity of A* and the corresponding 
nontrivial solutions of ()4.4|1 . 



5 Proofs of Theorems I2.2H2.4 



Now we are going to demonstrate the remaining part of our claims. 

Proof of Theorem \2/A We will give the proof for the case A* = A„(a_), the argu- 
ment for A=K = A„(a+) is similar. In accordance with the results of the previous 
section the eigenvalue A~(a, /), if it exists, it must be a root of the equation ()4.19j) . 
Let us prove first that there is a unique root which converges to A* as / — +00. 
Proposition 12.21 implies that the relation 

<j)+{xi, l,a^)= rc(A„ a_)e-2-^(A.)'e-'^f(^*)"i + 0{e-^P^^'^'), (5.1) 

holds in the norm of L2 (7a+ ) , hence by the definition of we have 

= rc(A„a_)e-2«f(^*)'(e-'^^(^*)"\0) + 0(e-2''(^*)'). (5.2) 

This formula in combination with the estimate ()4.12|) lead to the relation 

A22(A,0 = O(e-''^^(^*^'). (5.3) 

Since Tg is holomorphic w.r.t. A and has a small norm for large /, we infer that 
the left-hand side of the last equation is holomorphic in A. For a small 5 take the 
circle of those A such that |A — A*| = 5. In view of ()5.3|) the function A22 satisfies 
the estimate 1^4221 < 5 if / is large enough and |A — A*| = 5; by Rouche theorem 
it implies that the function A 1— > A — A* + v422(A, /) has the same number of zeros 
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in the disk {A : |A — A*| < 5} as the function A h-* A — A* does. The number S 
is arbitrary, so we can conclude that there is a unique root of the equation ()4.19p 
converging to A,,, as / — +00. As a consequence, there exists a unique eigenvalue 
of the operator H converging to A=k as / — > +00; we will denote this eigenvalue as 
X~{l,a). The estimate implies at the same time that 



\-{l,a) - A* = C(e 



-2K+(A*)^^ 



(5.4) 



Let us derive the asymptotic expansion ()2.8|1 for the eigenvalue A~(/, a). In 
order to do it, we will need to know the asymptotic behavior for A22 in a way 
more precise than ()5.Hj) . For the sake of brevity we will write shortly A instead of 
X~{l,a). The relations ()5.2|) together with the estimates ()4.12|) . ()5.4|) imply that 



A22{X, I) = T^{1 + Tg{X, I, a))- V;(-, 
= T^^:i;l)-T^T,{XJ,a)cl,:{;l) + Oi\\nfU 



-T^T,{X,,l,a)cj,;{;l) + 0[\X-X 



dX 



(5.5) 



-T4-T8-(A,/,a)0+(-,/,a_) + O (\X - X,\l\-<^^'^' + le'^^t^^*)' 



Taking into account the estimate ()3.14|) for ||Tg || and the relation ()5.1|) . we can 
proceed with the calculations obtaining 



A22{\1) 



TC 



(A, a_)e-2'^^("*)' (Tg- (A„ /, a+, a_)e-''^(^*)^S V'*) 



i2(7a_) 



+ O (|A - A.l/^e-^'^f^^*)' + e-2'(-f(^*)+p(^-))) 



(5.6) 



where we have denoted ip^:{x) = ipn{,x,a^). In view of the relation ()H.18j) the 
function Tg~(A,,, /, a+, a_)e'^i '•'^*-*^^ coincides with T-^ (X,l, a^, a-)u, where u is the 
solution to the problem with a = a^, X = A*, and / = e~'^^ (^O^i, it is clear 
that u{x) = t/(— xi, X2, A^K, a+), and in view of ()2.3j) . ()2.5j) . ()2.6j) we obtain 



(Tg (A^,, /,a+, a_)e 



L2ha_) 



c(A„a+)e-2''^(^-)'(e'^^(^*)^\^,) 



+ C(e 



= 7rc(A„ a+)c(A„ a_)/€+(A,)e-2<(^')' + 0(e"2p(^*)')- 
Substituting these identities into ()5.6|) . we finally arrive at the following formula, 

A22(A,/) = /.-(/, a)e-^«^(^*)' + O (|A - A,|/V^''^(^*)' + e-2'(''^(^*)+'^(^*))) , 

where /i~(/, a) is defined by ()2.9|1 . It allows us to rewrite the equation ()4.19|1 as 

(A - A,) (1 + 0(/2e-4«f(^*)')) = ;,-(/, a)e-^'^^("*)' + O {e~^ii4{^*)+pi>^*))' 
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expressing (A— A*) from here we get the asymptotic expansion ()2.8|) and the formula 

(E3D. 

Next we have to prove the asymptotic expansion for the eigenfunction asso- 
ciated with A~. The nontrivial solution of the equation ()4.4jl is given by ()4.14|) 
with c+ = and c_ = 1, i.e. as / = We substitute it into the relation 

n+ = T^^X^ , a+) and take into account ()5.2|1 . ()4.1(jj) : this yields 

= T3(A+,a+)/+ = O(IICII) = Oie-'^ti^'^^ 

which holds true in W2(Jla+) and in W2{S) for each S G Hq^. If A~ 7^ A=k, we 
obtain similarly with the help of Lemma 13.31 

Due to the equation ()4.19|) it follows that 

u_ = + C)(e-2«^(A.K) 

holds in W^2^(na) and W2{S) for each S G If A^ = A*, the last relation holds 
again; in order to prove it, it is sufficient to employ the identity 

M_ = r5(A*, a_)/_ - c_^* = r5(A*, a_)/_ - i^^. 

The relations obtained in this way together with ()4.3j) lead to ()2.1()|1 . □ 

Proof of Theorem \2.^ The general lines of the proof are similar to those of the 
previous one. According to the results of the previous section the eigenvalues of 
H converging to A=k are roots of the equation ()4.22j) . First we will check that the 
function at the left-hand side of this equation has either two simple zeroes or one 
second-order zero converging to A,,, as / — >• +00. 

To this aim, we need to estimate the functions Aij. Proposition 12.21 implies 

0- (xi, /, a+) = (-l)"^rc(A„ a+)e-2«^(^* V^^(^*)^^ + 0(e-2p(^*)'), (5-8) 
This formula together with allow us to conclude that 

A,,(A,0 = C?(e-=^"^('*^'), (5.9) 

hence for any small 6 we have the inequality 

|(A-A,)trA(A,/) + detA(A,/)| as |A-A*|=(5, 

if / is large enough. Since the functions A^j are holomorphic, by Rouche theorem 
this inequality implies that the function A 1—* D{\, I) := det ((A — A*)E + A(A, /)) 
has the same number of zeroes (with the order taken into account) as the function 
A I— >• (A — A*)^ does. The last function has A* as a second-order zero, of course, so 



A — A=K 



+ 0( 



(5.7) 
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it follows that the function D{-,1) has either two simple zeroes or a second-order 
zero, converging to A,,, as / — +00. In what follows we denote these roots as A^, 
the case of the second-order zero corresponds to the equality A"*" = A~. 

As it was established in the previous section, the nontrivial solutions of the 
equation ()4.4j) associated with the roots of ()4.22|1 are given by ()4.14|1 with the 
coefficients c± solving the system of linear equations ()4.2()j) . If the numbers A^ 
solve ()4.2H) . the system ()4.20|) has at least one nontrivial solution corresponding 
to A"^ and A~. 

Suppose that A^ 7^ A~. Then A^ are simple zeroes of the function D{-,1), 
and in view of the above discussion the system ()4.2()j) has exactly one non-trivial 
solution for A = A+ and A = A^. Hence in the case A'^ 7^ A^ the operator H has 
exactly two simple eigenvalues converging to A* as / — +00. 

Let us check that if the system ()4.20j) has two linear independent solutions 
referring to A = A^ it follows that A^ is a second-order zero of the function D{-,1). 
Indeed, two linear independent solutions exist if and only if 

An(A±, /) = A22{X^, /) = A, - A±, Ai2(A±, /) = A2i(A±, /) = 0. (5.10) 

The derivative of D{XJ) with respect to A equals 

|^(A, /) =2(A - A,) + (An(A, /) + A22(A, /)) 

+ (A-A.) (^(A,0 + ^(A,/) 
+ AuiX, 0^(A, - Ai2(A, 0^(A, 
+ ^2(A, 0^(A, - A2i(A, 0^(A, /). 

Substituting from ()5.10|) into this expression, we see that 

BD 

— (A, = as A = A±, 

thus A^ is a second-order zero. 

It is more complicated to check existence of a double eigenvalue of the operator 
i7 if A"*" = A~ =: A. It is equivalent to the fact that for A = A the system ()4.20|) 
has two linear independent solutions, and this in turn is equivalent to the relations 
()5.10|1 . Let us prove that they hold. Consider the boundary value problem 

(A + A)m = 0, X G n„ \ (7+ U 7_), u = 0, xedU, 



du 



dx-2 



du 

X2=+0 0X2 



(5.11) 



Here g± G -Z^2(7±) are arbitrary functions, and the parameter A is supposed to range 
in a small neighborhood of A* without coinciding with A* and A. This problem 
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is uniquely solvable provided we seek a L2(n)-solution to ()5.1H) . In a complete 
analogy with the proof of Lemma [4. II one can check easily that the problem (jS.lip 
is equivalent to the equation 

f + Ts{XJ,a)f = g, (5.12) 

where g = {g+,g-) G i^2(7a_) © -^^2(70+), while the solution u of ()5.11|1 is given by 

m(x, a, /) = M+(xi - l,X2, A, /) + M_(xi + /,X2, A, /), u± := T3(A, a±)f±. 

We can solve the equation ()5.12p in the same way as the equation fl4.1U|) . obtaining 
as a result that 

/ + -^$,+ + -^$; = G, G:={l + ni\,l,a))-'g. (5.13) 

Hence the function / is of the form 

/ = C+$+ + C^$; + G, (5.14) 

where C± = C-|-(A, /) are constants to be found. Denoting C := (^C-^ substi- 
tuting ()5.14|) into ()5.13|) . we obtain an equation for C, 

((A-A,)E + A(A,/))C = /i, h := (Zt'-g) ■ ^^-^^^ 
The solution of this system is given by Cramer's formula, 

Au{X, l)T^G - (A - A, + A22(A, /))T+G 



C+{\1) 



D{X,l) 

A2i(A, /)r+G - (A - A, + An(A, l))T^G 



(5.16) 



(5.17) 



DiXJ) 

Using now ()5.15j) and Lemma f3.3| we infer that 

u+{-, A, /) = - C+(A, /)^„(-, a+) + C+(A, /)T5(A, a+)<l>++ 
+ C_ (A, /)r5(A, a+)$;+ + r5(A, a+)G+, 

M_(-, A, /) = - C_(A, /)7/.„(-, a_) + C_(A, l)n{X, a_)$+_ 
+ C_(A, /)T5(A, a_)$;_ + T5(A, a_)G_, 

where $^_|_ and G± are the components of the vectors $f and G, 

$f = ($^+,$^_), G=(G+,G_). 

Since the number A is a second-order zero of D(-, /), we conclude from ()5.16|) that 
the coefficients C± have, in general, a second-order pole at A, and the same is true 
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for u±. Taking into account ()5.17|) . we conclude that the solution of ()5.11|) can be 
represented as 

u{x,X,l) =u+^{x,X,l)C+{X,l) +uZi{x,X,l)C4X,l) +0{1), X^X. (5.18) 

In a complete analogy with the proof of Lemma 13.31 one can check easily that 
the solution of the problem ()5.11|) has a simple pole at A. Hence the function 
u1li{x, X,l)C+{X,l) + uZi{x, X,l)C^{X,l) has a simple pole at A. For x from a 
neighborhood of 7+ this function satisfies due to ()5.16|) . ()5.17p the relation 

D{X, I) a, l)C+{X, I) + uZiix, X, /)C_(A, /)) 

= ((A - A, + A,2{X, l))T+G - A,,{X, OT^'G) - /, x^) + Oie-'^^^^*^'). 

Since A is by assumption a second-order zero of D{-,1), the obtained identity yields 
that 

A-A, + A22(A,/) = A2(A,/) = 0. 

Observing the behavior of the function u for x in the vicinity of 7_, one can prove 
in the same way that 

A-A, + An(A,/) = A2i(A,/) = 0. 

This completes the check of the relations ()5.10|) for A+ = A^ showing that in this 
case the operator H has a double eigenvalue converging to A=k as / ^ +00. 

We proceed to calculation of the asymptotic expansions for the root(s) of the 
equation ()4.22|) . Substituting the estimates ()5.9p into ()4.22|) we obtain 

A-A* = o(e-"^(^*)^). (5.19) 

This relation in combination with ()5.H1 . ()5.5|1 and the estimate ()3.14|1 imply that 

A22(A,/) = 0(/e-^"^(^*)'). (5.20) 

It is easy to establish an expression for An similar to ()5.5j) . which together with 
dEHI) and (fOTHl yield 

An(A,/) = 0(/e-^"^(^*)'). (5.21) 
Proceeding in the same way as in ()5.5|) we obtain a chain of relations, 

Au{x, I) = T+(i + T,{x, I, a))- V;(-, = ^M"!", + oiwnw iicii) 

= (0+(-, /, a_), ^^JU■, a+))L.(,.^) + 0{le-'^ti^*)iy 
Due to ()5.H) and ()2.3|) we have 

(0n(-,^,«-)>V'm(-,a+))L2(7a+) 
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where a) is given by ()2.12j) . Consequently, 

Ai2(A, /) = a)e-'"^(^*)' + 0(e-2''(^*)' + le'^'^ti^^^)^ (5.22) 
and in the same way one can show that 

A2i(A, /) = a)e-2''^(^*)' + C(e-2^(^*)' + le-^'^ti^^^), (5.23) 
The equation ()4.22p is equivalent to the following pair of the equations, 



2 

If c(A=K, a_)c(A*, a+) = 0, these equations together with ()5.20p - ()5.23p imply that 

A-A. = 0(e-2''(^*)' + /e-^"^(^*)')• 
which proves the asymptotic expansion ()2.1ip in the case a) = 0. 

Suppose on the contrary that c(A*, a_)c(A*, 0+) 7^ 0. In this case the function 
(All — A22Y + 4^12^21 is non-zero as A = A=k, and therefore its square root is 
holomorphic w.r.t. A. Using this fact and the relations ()5.2()|1 - ()5.23|1 . one can 
show easily in analogy with the similar argument for the equation ()4.22|) that each 
of the equations ()5.24|) has a unique root converging to A* as / ^ +00. Hence 
one of the roots of ()4.22|) satisfies the first of the equations ()5.24|) . while the other 
satisfies the other one. Substituting now from ()5.2()j) - ()5.23j) into ()5.24|1 . we arrive 
immediately at the asymptotics ()2.1H) . ()2.12|1 in the case a) 7^ 0. □ 

Proof of Theorem \2.4\ Let c be a nontrivial solution to the system ()4.20|1 . where A 
is A"*" or A~. Without loss of generality we may assume that \\c\\^2 = 1. Modifying 
()4.14|1 . we choose the corresponding nontrivial solution of the equation ()4.4j) as 
/ = —c^^^ — c^$~. In analogy with ()5.7|) we then obtain 

A — 

A — A* 

which holds true in W2(Jla-) and W^2^(S') for each S G Employing now the 

system ()4.20p we can write 

c+A2i(A, /) + c_A22(A, /) = -c_ (A - A,), 

hence 
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and in the same way one can prove that 

M+ = c+^^(-,a+) + 0(e-2''^(^*)'). 

in the norm of W2(Jla-) and 1^2^(5') for each S G The last two relations prove 
the sought formula ()2.15p . 

Suppose that A"^ = A", then ()4.20p has two nontrivial solutions, which means 
that (A — A*)E + A(A, /) = 0; we can choose these solutions as (c+,c_) = (—1,0) 
and (c+,c_) = (0, —1). Substituting these values into ()2.15|1 . we arrive at ()2.14|1 . 

Suppose that a) ^ 0. In view of ()2.1H) it implies that A"'"(A, a) 7^ A^(A, a), 
i.e. that A^(A,a) are simple eigenvalues. In this case the relations ()2.1H) and 
(jKTC|l . (jH:^ yield 

A± - A. + An(A±, = ±|^(/, a)|e-2-^(A0^(i + ^ q, 

, , (5.25) 

Ai2(A±, /) = /i(/, a)e-2''^(^*)'(l + 0(e-2«^(^')')) ^ 0. 

Since the matrix (A^ — A*)E + A(A^,/) has rank one, we can choose nontrivial 
solutions of ()4.20|1 as 

^± ._ V2(A^-A, + An(A^,/)) 



c 



c: 



v/(A±-A, + An(A±,/)2 + A?2(A±,0 
± ._ ^ V2Ai2(A^,/) 



v/(A±-A, + An(A±,/)2 + A22(A±,/) ■ 
In view of to ()5.25|) we then have 

= 1 + O(e-2'^^(^*)0, = Tsgn/i(/, a) + ^(e-^"^^^*)')- 
Substituting from here into ()2.15|) we arrive immediately at ()2.13p . □ 
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